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The spin transition in LaCoOa has been investigated within the density-functional theory -|- 
dynamical mean-field theory formalism using continuous time quantum Monte Carlo. Calculations 
on the experimental rhombohedral atomic structure with two Co sites per unit cell show that an 
independent treatment of the Co atoms results in a ground state with strong charge fluctuations 
induced by electronic correlations. Each atom shows a contribution from either a d® or a dJ state in 
addition to the main d® state. These states play a relevant role in the spin transition which can be 
understood as a low spin-high spin (LS-HS) transition with significant contributions 10%) to the 
LS and HS states of d® and dU states respectively. A thermodynamic analysis reveals a significant 
kinetic energy gain through introduction of charge fluctuations, which in addition to the potential 
energy reduction lowers the total energy of the system. 

PACS numbers: 71.27.-|-a, 71.30.-|-h, 71.10.Fd, 71.15.Mb 


The interpretation of the remarkably low temperature 
(~ lOOK) spin transition in LaCoOa (LCO) has been an 
intriguing research topic for decades (see Ref. [H for a 
review). A long standing debate regarding the origin of 
the transition is still ongoing with low-spin (LS, ^ 2969 ) 
to high spin (HS, ^ 2969 ) low-spin to intermediate 
spin (IS, tig Cg) models with concomitant orbital order¬ 
ing competing ii- The interest in this system has 
been recently boosted by the potential applications of 
this material and its doped phases in various optimized, 
environmental-friendly energy production domains Q. 
From the theoretical point LCO addresses one of the most 
challenging questions in solid state theory: the treatment 
of strongly correlated materials in which HS (IS) with¬ 
out long range magnetic order is stabilized [I| . Although 
the addition of static local correlations to the density- 
functional theory (DFT-I-U) improved the description 
of correlated materials with long range magnetic order, 
the more advanced formalism combining DFT with dy¬ 
namical mean-field theory (DMFT) @ has to be used to 
accurately treat systems in aj)aramagnetic state with lo¬ 
cal moments, like LCO [3, 0- Many-body calculations 
on LCO have been performed using DFT-f-DMFT Q 
and the variational cluster approximation [I 3 . In the 
first case, the contribution from LS, HS and IS was stud¬ 
ied and the spin state transition described as a smooth 
crossover from the homogeneous LS state into a non- 
homogeneous mixture of all three spin states. On the 
other hand, the VCA calculations showed that only the 
LS and HS states have appreciable weight in the density 
matrix over a wide temperature range The most 

recent attempts to settle the spin transition in LCO by 
means of the DFT-t-DMFT approach [3, S| used the nu¬ 


merically exact continuous-time quantum Monte Carlo 
(CT-QMC) [ll|. This methodology allows the inclusion 
of local dynamical effects and temperature, which is not 
possible within the inherently zero temperature DFT and 
DFT-bU. In both cases, a LS-HS transition scenario re¬ 
sults. No evidence of IS configuration was obtained, even 
upon going beyond the dP ionic picture Thus, in¬ 
cluding dJ and d® charged states to describe the local 
dynamics of the system, has allowed to interpret the spin 
state transitions as an LS (with few HS ions) to LS-HS 
short range ordered phase, with a subsequent melting of 
the order at higher temperatures. At room temperature, 
when 50% LS-HS population is expected, a Co(LS)-0- 
Co(HS) arrangement is anticipated. Most of the ab initio 
many-body calculations discussed so far have been per¬ 
formed assuming equivalent Co atoms, with exception of 
a couple of DFT-bU studies [1^ [l3|. In a groundbreak¬ 
ing study on a two band Hubbard model within DMFT 
Kunes and Kfapek 1^ showed that charge imbalance be¬ 
tween sites can occur on purely electronic grounds and 
possibly consititues an important piece in the LCO spin 
state transition puzzle. 


In this letter, we bolster this proposal by reporting on 
the first ab initio study of correlation driven charge and 
spin fluctuations in LCO within the DFT-I-DMFT for¬ 
malism for inequivalent Co atoms. Our results show that 
symmetry breaking of the Co sites creates a correlation- 
driven charge and spin fluctuation of purely electronic 
origin. 

LaCoOa is a perovskite system showing a small tetrag¬ 
onal distortion of the CoOe octahedra that varies with 
temperature, as determined from neutron diffraction [l5| . 
To describe it we use a fully ab initio approach start- 
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iiig from density-functional calculations of LCO using the 
experimentally determined structures for 5K, 300K and 
650K. We use a rhombohedral unit cell containing two 
formula units. The VAST code M with projector aug¬ 
mented wave basis sets (PAW) [171, and the Perdew- 
Burke-Ernzerhof [13 functional was employed for DPT 
calculations. A Wannier type construction using pro¬ 
jected local orbitals as described in [HI was applied 
to construct a local low energy model, which contains the 
on-site energies (crystal field splitting) and hoppings ex¬ 
tracted from DPT. We focus here on the 3d orbitals of Co 
exclusively, avoiding the possible ambiguities that arise 
from the use oi d + p orbital models in DPT-I-DMPT, 
see e.g. [H, [ 23 . Since the system is not perfectly cubic 
a straightforward Wannier construction produces a basis 
that retains some on-site mixing between the Co t 2 g and 
Cg orbitals. This local ^ 29-63 hybridization is mostly a 
consequence of this specific choice of orbital representa¬ 
tion, therefore we have performed a unitary transforma¬ 
tion after the Wannier projection to minimize it. The 
usual choices here are a rotation into the so-called “crys¬ 
tal field basis” or into a basis that renders the DPT oc¬ 
cupancy matrix Pp = (cjc^) diagonal on each atom, see 
e.g. Refs, li^ and (HI. The data presented were obtained 
using latter approach. We have, however, explored both 
approaches and found that they yield similar off-diagonal 
elements in the frequency dependent DPT Green’s func¬ 
tion G^(iuJn) and the choice has no qualitative effect on 
any conclusions drawn in what follows. 

On the level of DPT we find a set of three orbitals 
very close in energy (two and one aig, that we will for 
brevity call t 2 g) and two orbitals (Cg) higher in energy 
by about A ~ 1.5-1.65eV. We solve an effective multior¬ 
bital Hubbard model within DMPT using a hybridization 
expansion CT-QMC solver using density-density interac¬ 
tion terms. To treat the two Co atoms in the unit cell 
independently we employ the so-called inhomogeneous or 
real space DMPT [2^ , where we have to solve an Ander¬ 
son Impurity model for each correlated atom a in the unit 
cell and the atoms effectively interact via their respective 
baths of conduction electrons. The generalised five band 
Hubbard model including the Coulomb interaction and 
a double counting correction then contains the following 
terms 


H — +-'S^U'^-n h / 111 

~ 2 / ^ aij ' ' ^ocj,<7' •> \^) 

a,a' 

where is the (Kohn-Sham) DPT Hamiltonian, fiai^cr = 
and is the creation operator of an elec¬ 
tron on site a in Wannier state i and spin a. The dou¬ 
ble counting, oc amounts to a shift of the chemi¬ 

cal potential for the Co 3d shell and is determined self 
consistently. To obtain the Coulomb interaction matrix 
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FIG. 1. (color online) Phase diagram as a function of the 
lattice structure temperature and the Hunds rule coupling J 
for U = 3eV. The colored (shaded) parts illustrate different 
regions in the phase diagram, white indicating LS and blue 
(gray) HS. (a) shows the homogeneous phase exhibiting the 
LS, MS and HS states as indicated by the data points. MS 
refers to a mixture of LS and HS states, (b) shows the data 
points and possible phase regions for the calculation with in¬ 
equivalent Co atoms. The colored (shaded) areas are again 
only meant as an illustration, the checkerboard pattern now 
indicating the LS-HS phase where strong charge fluctuations 
are present. 


we employ the parametrization via the Slater integrals 


27| connected to the average direct and exchange cou¬ 
plings U and J: F° = U, J = 1 / 14 (F 2 + and 
P"* = 0.625P^. The aforementioned basis transforma¬ 
tion was applied here as well. 


One important element to properly describe the sys¬ 
tem within DPT-I-DMPT is the appropriate choice of the 
Coulomb interaction. In the case of LCO, the problem is 
more delicate since after the LDA-I-U calculations that in¬ 
troduced the LS-IS model [3 , it has been widely believed 
that it is a strongly correlated electron system with an 
on-site Coulomb interaction oi U = 8eV. Recently, con¬ 
strained DPT results proposed a value oi U = 6eV for 
the Co 3d shell in a d -b p model [3 • Such high values 
for the Coulomb interaction might be appropriate for a 
d + p model, are however too large for d only calcula¬ 
tions. Since we are not aware of any ab initio estimates 
for the Co 3d shell only, measurements of the excitation 
gap in LCO have been used as a guide. A gap size of 
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TABLE I. Most probable many-body configurations for the 
300K structure with two inequivalent Co atoms as a function 
of the Hunds coupling J obtained from the analysis of the 
CT-QMC imaginary time evolution (in %). Numbers missing 
to or sums exceeding 100% are due to minor contributions of 
other atomic states and rounding. 


about 0.9eV was measured from photoemission and ab¬ 
sorption spectra [1^, while smaller values, 0.6eV [2^ and 
0.3eV [1^ were obtained from optical measurements. We 
were able to obtain a charge gap ~ 1 eV with a value of 
U = 3eV, which will be used in all subsequent calcula¬ 
tions. 

Since the value of the Hunds rule coupling J is, for 
a fixed U and crystal field splitting A, the critical pa¬ 
rameter for the spin state transition 31| we have cal¬ 
culated the system at different values of J. To ac¬ 
count for the temperature, the calculations were per¬ 
formed for the experimental crystal lattice structures de¬ 
termined at the temperatures 5K, 300K and 650K [l5| . 
In the QMC solver we used the calculation temperatures 
116K, 290K, 580K (equivalent to the inverse tempera¬ 
tures P = 100,40, 20eV“^) for these structures respec¬ 
tively. 

In a first approximation, the two Co atoms were con¬ 
strained to be in the same charge and spin state. The 
calculations show that for the three crystal structures 
and their respective crystal fields a spin state transi¬ 
tion occurs at about the point where A ^ 2J, i.e. for 
J 0.75 — 0.8eV. In Fig. [T^ we have plotted our data 
for U = 3eV as a function of the ’’lattice temperature” 
and J. We can see that there is a crossover region, that 
we call mixed spin (MS), between the homogeneous LS 
(white region) and HS (blue region) phases. The tran¬ 
sition in the MS region is governed by an increased ad¬ 
mixture of the HS contribution to the LS state. The 
hybridization expansion CT-QMC solver allows for an 
analysis of the local eigenstates contributing to the par¬ 
tition function during the imaginary time evolution. The 
states observed here are almost pure LS or HS (> 90% 
probability), i.e. ^ 1=2 no contribution 

from any IS (dg^i) states. Small (^ 3%) contributions 
of found in the LS and HS states 

respectively. 

In a second step, the constraint of equivalent Co atoms 
in the unit cell was removed. The calculations show that 
within this assumption, strong charge fluctuations be- 
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FIG. 2. (color online) Orbitally resolved spectral functions 
for LCO for the 300K crystal structure at solver temperature 
of 290K (/? = 40eV“^). (a) The homogeneous HS (J = 0.9eV, 
top) and LS spectra (J = 0.6eV, bottom) with atomic states 
given as an inset, (b) Results in the asymmetric Co configu¬ 
ration for the values of J indicated. The largest and second 
largest many-body contributions in the LS-HS ordered phase 
with J = 0.75 eV are again given for a simplified octahedral 
crystal field. 


tween the two ions develop. This can happen sponta¬ 
neously via noise introduced by the QMC procedure, but 
we have also introduced a small difference in the levels 
(/iP®-/i^^ = 0.02eV) in the first DMFT iteration to ren¬ 
der the two atoms explicitly inequivalent 32| . More con¬ 
cretely, a contribution of dP and dJ states to the nominal 
d^ average charge develops as a function of J. Thus, at 
small J the two atoms both converge to the LS configura¬ 
tion as before, but with increasing J charge fluctuations 
between the two atoms occur, see the LS-HS region in 
Fig. [TJd. Regarding the spin configuration in the LS-HS 
phase one atom will be in a predominantly LS and the 
other in a predominantly HS state. This is accompanied 
by the occupancies of the respective Co 3d shells to devi- 
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ate from N^d = 6.0, see Tab. HI Consequently, the QMC 
partition function shows sizeable contributions of 
on atom 1 and ^ 5 ^ 3/2 atom 2, respectively. Other 
theoretical results and the interpretation of experimen¬ 
tal data im indicate that such a state exists in LCO 
at room temperature. From Fig. [T]one can see that the 
spin transition can be studied as a function of the Hunds 
coupling J and of the temperature. This implies, that 
the transition can be driven only by electronic means, as 
shown by previous model calculations 3^. In the follow¬ 
ing, the detailed electronic structure configuration will be 
investigated as a function of J assuming the 300K crystal 
structure, which exhibits all relevant features present in 
the whole temperature range. 

The evolution of the Co 3c? spectra as a function of J 
is given in Fig. [2] The orbitally resolved spectral func¬ 
tion (obtained via maximum entropy 3^ from the QMC 
Green’s function) is shown for the homogeneous case for 
the LS and HS states of the 300K crystal structure in 
Wa.. The strongest overall changes are visible for the un¬ 
occupied part of the spectra, suggesting that experimen¬ 
tal techniques addressing those states will be relevant 
to understand the system in more detail [s^. The low- 
spin state in Fig. [5^ is closest to the DFT solution, the 
strongest modification is the rigid upward shift of the eg 
bands and, as a consequence, the gap opening between 
t 2 g and Cg states. This is in accordance with combined 
DFT and cluster calculations as well as recent QMC 
Q. The formation of local moments in the higher tem¬ 
perature HS states leads to the appearance of additional 
features in the spectrum. As a result the gap changes 
its character from t 2 g — Cg to t 2 g — Qg with incoherent 
t 2 g excitations on both gap edges. The occupied parts 
of the spectra exhibit a transfer of spectral weight away 
from the strong t 2 g excitation peak towards higher bind¬ 
ing energies as the LS to HS crossover commences, see 
Ref. 0 and references therein. The spectral function 
for the asymmetric Co configuration is displayed in Fig. 
do for the LS-LS (J = 0.6eV), the LS-HS (J = 0.75eV) 
and the HS-HS (J = 0.9eV) arrangements. Again, as the 
transition from LS-LS to HS-HS commences the first t 2 g 
excitation peak is largest for the LS-LS state, reduced in 
the LS-HS and even more so in the HS-HS state. Also, 
the progressive reduction of the gap expected from ex¬ 
perimental studies is reproduced 33|, IST], ■ 


The tendency of the system to introduce charge fluctu¬ 
ations can be understood by analyzing the total energies 
within DMFT. The kinetic energy is accessible within 
QMC as Akin = —T{k), where T is the temperature and 
{k) is the mean order of the QMC histogram, whereas 
the potential energy can be obtained from the double oc¬ 


cupancy as i ^ 


r/co- /" 

a,2J,CT,<7' ^ OilJ 


y>') l39|,l40|. As a 


consequence of the charge fluctuations the “critical” J 
for the transition on one atom is lowered (via the dP ad¬ 
mixture) allowing to introduce one HS site into the cell. 



FIG. 3. (color online) Analysis of the total energy contri¬ 
butions for the 300K structure, (a) Orbitally averaged dou¬ 
ble occupancies in the homogeneous and inhomogeneous so¬ 
lutions (empty and solid symbols respectively) along with the 
corresponding Coulomb interaction term. The potential en¬ 
ergy (obtained from the double occupancies) is lower in the 
LS-HS phase as compared to the homogeneous LS-LS phase, 
see text, (b) Histogram of the expansion order for 290K 
(/3 = 40eV“^) of the QMC diagrams contributing on one 
atom to the fermionic trace for the homogeneous (shaded) 
and the LS-HS phase (solid symbols). Its average is propor¬ 
tional to the kinetic energy. Numerical values correspond to 
the histograms shown. 


The HS atom is then in a state of predominantly 
character with a contribution of a ^ 5 ^ 3/2 state. The LS 
atom on the other hand is now in a with an admix¬ 
ture of which is higher in energy compared to the 

pure dg^Q LS state. Balancing the energy loss of the LS 
atom with the energy gain brought about by the HS atom 
the mixed LS-HS arrangement constitutes a net gain in 
potential energy. The total energies can be analyzed in 
detail using Fig. [3] The bulk of the potential energy gain 
can be understood as a reduction of the on-site double 
occupancies (ni,CTni,-cr), see terms oc in Fig. |3^. In 

the homogeneous LS state these terms are close to their 
maximum of 1 for the Qg and zero for the eg states, thus 
giving the average value of 0.6 per orbital, while dropping 
to about 0.2 per orbital in the HS state (empty symbols in 
Fig. |3K). In the LS-LS state both atoms contribute with 
~ 0.6 Uii per orbital, while in the LS-HS state one atom 
contributes ~ 0.6 Uu, the other only ~ 0.2 Ua (solid 
symbols in Fig. [3^). The net gain due to the contri¬ 
butions of the “anti-Hund” double occupancies oc 
and the energy loss due to the “Hund” terms oc al¬ 
most cancel. Using the proper numerical values for the 
double occupancies and interaction matrices we obtain 
the difference AApot = Ri —2.4eV as- 
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suming J = 0.7eV. At large enough J the Hunds rule 
energy overcomes the crystal field and the HS-HS state 
becomes accessible. 

The QMC histogram —Akin) shown in Fig. [Sb illus¬ 
trates beautifully the appearance of charge fluctuations 
and can be interpreted using a reasoning similar to Ref. 


4l|. The QMC histogram for the pure LS and HS states 


shows the usual one peak structure, HS being more itin¬ 
erant and exhibiting a shoulder in addition to the prin¬ 
cipal peak, while for the LS-HS state a clear two peak 
structure is observed. The first peak indicates standard 
processes that also occur for a single site, i.e. hopping 
from site a to site a, while the second peak is of non-local 
origin and a consequence of the introduction of a second 
site a'. These processes have a larger kinetic energy (by 
absolute value) and thus are more itinerant than the ho¬ 
mogeneous LS and HS phases. Concretely, both atoms 
in the LS-HS phase have a kinetic energy even larger (by 
absolute value) than in the HS homogeneous phase. This 
is a hallmark of the strong charge fluctuations between 
the two sites and constitues an additional gain in energy 
not accounted for in a single site picture. 

The proposed mechanism can easily explain the low 
transition temperature observed in LCO on purely elec¬ 
tronic grounds, since already for the 5K structure at 
solver temperature of 116K the LS-HS phase is firmly 
established. Nevertheless, the discrepancy between lat¬ 
tice temperature and solver temperature distorts the 
phase diagram. Besides, some of the approximations 
used might influence the results. In future work it would 
be interesting to investigate how additional terms in the 
Coulomb interaction (spin-flip and pair-hopping) influ¬ 
ence the observed behaviour. In addition, the inclusion 
of the O 2p states into the interacting model with a 
Coulomb interaction on the p states (App) and between 
the d and p states {Upd) would be desirable, but is at 
the moment too expensive and an approximation of the 
Coulomb interaction difficult. Another interesting point 
would be how the charge and spin fluctuations found here 
couple to phonons, see e.g. [^. Finally, a more complex 
order than the reported checkerboard phase might also 
build up. 

In summary, ab initio DFT-l-DMFT calculations for 
the two atomic unit cell of LaCoOa show that upon treat¬ 
ing the Co atoms independently strong charge fluctua¬ 
tions develop in the system. As a consequence, the spin 
transition can be understood as a transition from a homo¬ 
geneous LS to a LS-HS ordered state with strong charge 
fluctuations and subsequently into a homogeneous HS 
state. This provides a novel understanding of the system 
in which the charge fluctuations are present in the sys¬ 
tem from first principles. Angle resolved photoemission 
experiments as a function of the temperature and time- 
resolved XAS studies with femtosecond resolution should 
allow to verify the proposed model. 
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